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Abstract 

We present a study of three-particle correlations among a trigger particle and two associated 
particles in Au + Au collisions at ^Jsnn = 200 GeV using a multi-phase transport model (AMPT) 
with both partonic and hadronic interactions. We found that three-particle correlation densities 
in different angular directions with respect to the triggered particle ('center', 'cone', 'deflected', 
'near' and 'near-away') increase with the number of participants. The ratio of 'deflected' to 'cone' 
density approaches to 1.0 with the increasing of number of participants, which indicates that 
partonic Mach-like shock waves can be produced by strong parton cascades in central Au+Au 
collisions. 
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I. Introduction 

Ultra-relativistic heavy ion collisions may 
provide conditions sufficient for the for- 
mation of a deconfined plasma of quarks 
and gluons [1]. Experimental results from 
RHIC indicate that a strongly-interacting 
partonic matter (termed sQGP) has been 
created in the early stage of central Au 
+ Au collisions at y/s^N = 200 GeV at 
RHIC [2]. Jet-like azimuthal correlation 
is one of the important hard probes to 
explore the natures of the newly formed 
matter. The disappearance [3] and re- 
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appearance [4] of back-to-back high trans- 
verse momentum (pr) particles from jets 
have been proved to result from the in- 
teractions between jet-partons and the 
hot and dense medium created in central 
Au+Au collisions. Recently, an interesting 
splitting of the away side peak has been 
observed in the di-hadron azimuthal angle 
(A</>) correlation distribution between soft 
associated particles and high p T trigger 
particles in central Au + Au collisions at 
RHIC [5,6,7]. Such a double peak structure 
on the away-side is consistent with pref- 
erential conic emission of particles from 
jets and/or shock-wave induced collective 
motion from jet-medium interactions. We 
will refer to the observed double peaks on 
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the away-side as the Mach-like structure 
without necessary implication on the dy- 
namical mechanism. 

Several theoretical interpretations about 
the Mach-like structure have been pro- 
posed. For instances, Stocker et al. pro- 
posed the Mach-like structure from jets 
traversing the dense medium as a probe of 
the equation of state (EOS) and the speed 
of sound in the medium [8]. Casalderrey- 
Solana and Shuryak et al. argued a shock 
wave generation because jets travel faster 
than the sound in the medium [9]. They 
fitted the broad splitting structure on the 
away side in di-hadron azimuthal correla- 
tion with a Mach-cone shock wave mech- 
anism. Vitev has shown that the cancella- 
tion of collinear bremsstrahlung in QCD 
medium can lead to large angle emission of 
gluons [10]. Koch and Wang et al. used a 
Cherenkov radiation model with negative 
dispersion relation to produce the Mach- 
like structure [11]. In Ref. [12], Armesto 
proposed that the medium-induced gluon 
radiation could be affected by the collec- 
tive flow in the medium. It has also been 
argued by Miiller et al. that a Mach-like 
structure can appear via the excitation of 
collective plasmon waves by moving color 
charges associated with the leading jet [13]. 
Renk and Ruppert found that in order to 
reproduce the experimental data a large 
fraction (about 90%) of the lost energy of 
jet has to be channelled to excite a shock 
wave in a dense medium at a soft point of 
EOS [14]. Satarov et al. investigated Mach 
shocks induced by partonic jets in expand- 
ing quark-gluon plasma [15]. However, 
Chaudhuri and Heinz reported no observa- 
tion of Mach-like structures in di-hadron 
A(p correlations from jet quenching dynam- 
ically in a hydrodynamic QGP fluid [16]. A 
consistent dynamical picture for the gener- 
ation of the Mach-like structure in particle 
correlations has yet to emerge and further 
investigations are needed. 

In order to shed light on the puzzle of the 
dynamical origin of the splitting structure 
on the away-side, three-particle correla- 
tion has been proposed to look at the 
multi-particle correlation in the emission 
pattern of particles. The di-hadron corre- 



lation cannot distinguish different emis- 
sion scenarios since correlation only deals 
between emitted and the trigger particle. 
However the three-particle correlation is 
capable of distinguishing the different sce- 
narios when simultaneous emission of two 
particles are investigated with the trigger 
particle. If the splitting structure of away- 
side is from large angle gluon emission or 
deflection due to strong collective flow in 
an event, the two associated particles will 
be clustered in a narrow cone on a single- 
side of the away-jet direction. However, if 
the production mechanism is Mach-cone 
shock wave or Cherenkov gluon radiation, 
the partons in the shock-wave front or 
Cherenkov gluons will be emitted conically 
around the away-side jet center in single 
event. In this case, the two associated par- 
ticles can be simultaneously on both sides 
of the A<p = 4>assoc — 4>tri g distribution with 
respect to the opposite direction of the 
trigger particle. Experimental studies of 
the three-particle correlations have been 
reported by both the STAR [17,18,19] and 
the PHENIX [20] collaborators. 

In our previous work, we reported obser- 
vation of Mach-like structure in di-hadron 
correlations from Au+Au collisions using 
a multi-phase transport model (AMPT) 
where both partonic and hadronic in- 
teractions are included [21]. Both par- 
ton cascades and hadronic rescatterings 
can produce apparent di-hadron correla- 
tions with Mach-like structures. But the 
hadronic rescattering mechanism alone 
cannot reproduce the observed experimen- 
tal amplitude of Mach-like structure on 
the away-side, which indicates that parton 
cascade processes are indispensable. How- 
ever, detailed dynamical mechanisms for 
the Mach-like structure still await to be 
determined. In this Letter, we present a 
study of three-particle correlation among 
one trigger particle and two associated par- 
ticles in Au + Au collisions at ^/sNN = 200 
GeV with the AMPT model. Three par- 
ticle correlations in regions of azimuthal 
angular directions of 'cone', 'deflected', 
'center', 'near' and 'near-away', which will 
be defined later, will be presented for Au 
+ Au collisions from AMPT. With de- 



2 



creasing number of participants, 'center' 
correlations become more dominant, and 
'cone' and 'deflected' correlations seem to 
disappear. Our results indicate that the 
three-particle correlations in central col- 
lisions are mainly produced by partonic 
Mach-like shock wave effect, while in pe- 
ripheral collisions deflected jet effect also 
contributes to the Mach-like structure. Ef- 
fects of hadronic rescatterings and parton 
cascades on three-particle correlation are 
also investigated. 

II. Brief Description of the AMPT Model 
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AMPT model [22] is a hybrid model which 
consists of four main processes: the initial 
conditions, partonic interactions, the con- 
version from partonic matter into hadronic 
matter and hadronic interactions. The ini- 
tial conditions, which include the spatial 
and momentum distributions of minijet 
partons and soft string excitations, are ob- 
tained from the HIJING model [23]. The 
excitation of strings will melt strings into 
partons. Scatterings among partons are 
modelled by Zhang's parton cascade model 
(ZPC) [24], which at present includes only 
two-body scatterings with cross section 
obtained from pQCD calculation with 
screening mass. In the default version of 
AMPT model (we briefly call it as "the 
default AMPT" model) [25], partons are 
recombined with their parent strings when 
they stop interactions, and the resulting 
strings are converted to hadrons using the 
Lund string fragmentation model [26]. In 
the string melting version of the AMPT 
model (we briefly call it as "the melt- 
ing AMPT" model) [27], a quark coales- 
cence model is used to combine partons to 
form hadrons. Dynamics of the subsequent 
hadronic matter is then described by A 
Relativistic Transport (ART) model [28]. 
Details of the AMPT model can be found 
in a recent review [22]. Previous studies 
[22,27,29] demonstrated that the partonic 
effect cannot be neglected and the melting 
AMPT model is much more appropriate 
than the default AMPT model in describ- 
ing nucleus-nucleus collisions at RHIC. 
In the present work, the parton interac- 
tion cross section in the AMPT model is 
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Fig. 1. Three-particle correlations in the top 
10% central Au+Au collisions at ^/snn = 
200 GeV from the melting AMPT model 
with hadronic rescattering. (a): Raw signal, 
(b): Hard-soft background, (c): Soft-soft back- 
ground, (d): Random background. 

assumed to be 10 mb consistent with pre- 
vious calculations [22,29]. 

III. Analysis Method 

The mixing-event technique has been used 
in our three-particle correlation analysis. 
The p T window cuts for trigger and as- 
sociated particles were selected as 2.5 < 
p% ig < 4 GeV/c and 1.0 < p a T ssoc < 2.5 
GeV/ c, respectively. Both trigger and asso- 
ciated particles were required to be within a 
pseudo-rapidity window of | ?7 1 < 1.0, where 
r) is the pseudo-rapidity of hadrons in the 
center-of-mass frame of Au+Au collisions. 
In the same events, raw 3-particle corre- 
lation signals in A0i — (f>i — (fttrig versus 
A</>2 = 4>2 — 4>tHg were histogrammed. Fig- 
urel(a) shows the raw 3-particle correla- 
tion distribution in the top 10% central 
Au+Au collisions at ^snn = 200 GeV in 
the melting AMPT model with hadronic 
rescattering. Three classes of background 
contributions are expected to contribute to 
the raw signal. The first one is the hard- 
soft background which comes from a jet- 
induced trigger-associated pair combined 
with a background associated particle from 
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Fig. 2. Background subtracted 3-particle cor- 
relations in the top 10% central Au+Au col- 
lisions at ^/snn = 200 GeV for the melt- 
ing AMPT model with hadronic rescattering. 
(a) and (b) : Background subtracted 3-particle 
correlations (—1 < A<j>i t 2 < 5.28). (c) and (d): 
background subtracted 3-particle correlations 
(1 < A<pi t 2 < 5.28) from the selected regions 
(1,2,3,4,5). The azimuthal angular regions are 
defined in panel (a) - (1): 'center' region; (2) 
and (3): 'cone' regions; (4) and (5): 'deflected' 
regions; (6): 'near' region; (7) and (8): 'near- 
away' regions. 

bulk medium. We reproduced it by mix- 
ing a trigger-associated pair with another 
associated particle from a different event 
(Figure 1(b)). The second one is soft-soft 
background which comes from an associ- 
ated particle pair combined with an un- 
correlated trigger particle. We constructed 
this background by mixing an associated 
particle pair from one event with a trig- 
ger particle from a different event (Figure 
1(c)). The third one is a random combi- 
natorial background, which was produced 
by mixing a trigger particle and two as- 
sociated particles respectively from three 
different events (Figure 1(d)). We required 
that the mixed events are all from very 
close collision centralities which can be de- 
termined by impact parameters in simula- 
tions. In order to subtract the backgrounds 
from the raw signals, we set the signal at 
0.8 < |A0i )2 | < 1.2 to be zero. Figure 2 
(a) and (b) give background subtracted 3- 



particle correlations in the top 10% central 
Au+Au collisions at ^/snn = 200 GeV in 
the melting AMPT model which includes 
hadronic rescattering. In order to observe 
the 3-particle correlations among a trigger 
particle and two away-side associated par- 
ticles clearly, the 3-particle correlations in 
1 < A0 1;2 < 5.28 region are shown with an 
expanded scale in Figure 2 (c) and (d). 

IV. Results and Discussions 

We divide the three-particle correlation 
distribution into several regions based on 
the possible origin of the particle emission 
pattern as shown in Figure 2a. The first one 
is 'center' region (|A0i 2 — n\ < 0.5) where 
the three-particle correlation mainly comes 
from one trigger particle and two associ- 
ated particles in the center of away side. 
The 'center' correlations represent pene- 
tration ability of away-side jet. The second 
one is 'cone' region (|A0i — (tt ± 1)| < 0.5 
and \A(p2 — (tt =f 1 ) | < 0.5) where three- 
particle correlation would form splitting 
peaks in di-hadron A0 correlation due 
to a conical emission pattern from away- 
side jet. It was predicted that this conical 
emission may be produced by a Mach- 
cone shock wave effect when a jet propa- 
gates faster than the speed of sound in the 
medium creating shock wave front in the 
cone region. The third one is 'deflected' 
region (|A</>i j2 — (tt ± 1)| < 0.5) where as- 
sociated particles are emitted in the same 
side- ward region of the away-side jet in one 
event. The 'deflected' region three-particle 
correlations can also yield splitting peaks 
on the away-side of two-particle correlation 
distribution because though within one 
event the away-side jet is deflected to one 
side only, but inclusively with many events 
both sides of the jet direction can be pop- 
ulated. The fourth region is the 'near' area 
(|A0i j2 | < 0.5) where three-particle cor- 
relation represents the correlation among 
trigger particle and associated particles 
on near side of the trigger direction. The 
fifth one is 'near-away' correlation region 
(1 < A0i, 2 < 5.28 and jA0 2 ,i| < 0.5), 
which reflects the correlation among trig- 
ger particle, one associated particle on 
near side and another associated particle 
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on away side. The five regions have been 
marked with different numbers in panel 
(a) of figure 2 for clarity. We will examine 
three-particle correlations in the above five 
regions. 

Figure 3 shows three-particle correlation 
distribution in the (1 < A^i^ < 5.28) area 
from Au+Au collisions at ^/snn = 200 GeV 
with different centralities using the melting 
AMPT model, and p+p collisions at ^snn 
= 200 GeV using the default AMPT model 
before and after hadronic rescattering. 
Here we chose the default AMPT model 
to simulate p+p collisions at ^/snn — 200 
GeV, since the string melting mechanism 
has little effect on p+p collisions which has 
also been demonstrated previously [22]. 
Three-particle correlations in all 'center', 
'deflected' and 'cone' regions can be ob- 
served in central Au+Au collisions with 
the melting AMPT model regardless of 
the inclusion of hadronic rescatterings. As 
the collisions become more peripheral, the 
'deflected' and 'cone' region correlations 
gradually disappear until only the 'center' 
correlations remain in the most peripheral 
Au+Au collisions and p+p collisions. 

In order to quantitatively express three- 
particle correlation strength in these differ- 
ent regions, region-averaged three-particle 
correlation density p is defined according 
to the following equation: 
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The top panel of Figure 4 shows three- 
particle correlation densities p in different 
regions as a function of N part (number 
of participants) for Au+Au collisions at 
^s^n = 200 GeV in the melting AMPT 
model before and after hadronic rescatter- 
ing. 

Our results show that three-particle cor- 
relation densities decrease after hadronic 
rescattering process, which indicates 
hadronic rescatterings could weaken three- 
particle correlation strength. However di- 
hadron correlation is almost unchanged in 
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Fig. 3. Background subtracted seg- 
mental 3-particle correlation areas 
(1 < A0i2 < 5.28) in different centralities 
in Au+Au collisions at ^snn = 200 GeV in 
the melting AMPT model ((a)-(j)), as well as 
p+p collisions at ^/snn = 200 GeV in the 
default AMPT model ((k)-(l)). The left col- 
umn from (a) to (k) shows the results before 
hadronic rescattering (briefly named as "melt 
before" or "default before") and the right 
column from (b) to (1) shows the results after 
hadronic rescattering (briefly named as "melt 
after" or "default after"), (a) and (b): 0-10%; 
(c) and (d): 10-20%; (e) and (f): 20-40%; (g) 
and (h): 40-60%; (i) and (j): 60-80%. (k) and 
(1): p+p collisions. 
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Fig. 4. The correlation density analysis for 
Au+Au collisions at ^/snn = 200 GeV in 
the melting AMPT model before and after 
hadronic rescattering. Top panel: the aver- 
age three-particle correlation densities p at 
different regions as a function of N part ; 
Bottom panel: ratios of average three-par- 
ticle correlation density ('center '/'deflected' 
and 'center '/'cone') as a function of N part ; 
The insert of the bottom panel: ratio of av- 
erage three-particle correlation density ('de- 
flected'/'cone') as a function of N part . Note 
that some points have been shifted slightly in 
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this Pt window selection in our previous 
work [21]. Such a difference is indicative 
of enhanced sensitivity to hadronic rescat- 
terings in the three-particle correlations in 
comparison to the di-hadron correlations. 

The bottom panel of Figure 4 shows 
the centrality dependences of two ra- 
tios, namely the density ratios of 'cen- 
ter'/'deflected' and 'center' /'cone', in 
Au+Au collisions at ^/s^vv = 200 GeV. 
Both ratios fall from above 2.0 in periph- 
eral collisions to near 1.0 in central colli- 
sions with the increasing of N part , which 
indicates that the strengths of particle 



emission in the 'cone' and in the 'deflected' 
regions increase dramatically in central col- 
lisions relative to the particles in the 'cen- 
ter' region. Since the 'center' correlation 
reflects the ability of 'punch — through' for 
the backward jet, our results indicate that 
the backward jet can maintain the original 
jet direction well in peripheral collisions 
while in central collisions many particles 
are emitted in the 'cone' and the 'deflected' 
directions, away from the original jet di- 
rection. 

In the insert of Figure 4, the ratio of 
'deflected'/'cone' slightly decreases with 
Np ar t and approaches 1.0 in central col- 
lisions. The Mach-cone shock wave and 
the Chrenkov gluon radiation scenarios 
predicted almost equal strength in the 
three-particle correlations in the 'deflected' 
(7r±D,7r±D) and 'cone' (7r±D,7r+D) re- 
gions, where D is the splitting param- 
eter of away side (i.e. half distance be- 
tween two peaks on away side in di-hadron 
A(p correlation function). Our observed 
three-particle correlations in the central 
Au+Au collisions from the AMPT model 
are consistent with these model predic- 
tions. Such a consistency may be related 
to the hydrodynamic-like behavior in the 
AMPT model due to strong parton-parton 
couplings and interactions. 

More comments on the origin of the three- 
particle correlations in the AMPT model 
are in order. The melting AMPT model 
was shown to produce good descriptions of 
elliptic flow of identified hadrons and even 
yielded the correct mass ordering of elliptic 
flow [27,29], which has been considered an 
important feature of hydrodynamics mod- 
els. Such an agreement can be attributed to 
the large parton-parton interaction cross 
section in the AMPT model, which leads 
to strong parton cascades that couples 
partons together inducing the onset of hy- 
drodynamical behavior [30]. However, in 
another hydrodynamic model [16] the sig- 
nal of Mach-cone shock waves can hardly 
be observed in the di-hadron correlations. 
It appears that the large strength of par- 
ton cascades and coupling of partons as 
described in the AMPT model bring about 
the conic emission pattern on the away-side 



6 



melt after(0-10%) 



S 1.6 




1.5 2 2.5 3 3.5 4 4.5 5 

Aflt(rad) 



1 1.5 2 2.5 3 3.5 4 4.5 5 



Fig. 5. Background subtracted three-particle 
correlations in selected (1 < A<pi t 2 < 5.28) re- 
gions ((a) and (b)) and away-side di-hadron 
correlations ((c) and (d)) in the top 10% 
Au+Au collisions at y / s/v r /v = 200 GeV in 
the melting AMPT model (left column) and 
the default AMPT model (right column) after 
hadronic rescattering. 

prominently. The linearized hydrodynami- 
cal approximation may not be adequate for 
the strong jet-medium interaction region 
where the medium also experiences rapid 
variation of energy density and without 
sufficient thermalization [9]. On the other 
hand, the observed three-particle correla- 
tions may partly stem from deflected jets 
(represented by pde ^ ected — 1 ) in periph- 
eral collisions where 'center' correlation 
becomes dominated. In the AMPT model 
there is no inclusion of large angle gluon 
bremsstrahlung mechanism [10] which may 
also play a role in real collisions. In addi- 
tion, we note that the backward jet may 
also be distributed over a wide rapidity 
range [10,31] beyond our narrow r] win- 
dow cut. In our model, we used LO pQCD 
cross sections from HIJING model for the 
minijet production, which has successfully 
described the suppression of back-to-back 
jets [32]. Our selection of the rj window cut 
was to match the detector acceptance of 
the RHIC experiments. 

In addition, we studied the effect of parton 
cascades on three-particle and di-hadron 
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Fig. 6. The px dependences of nuclear mod- 
ification factor R cp of charge hadrons for 



0-10%/40-60% in Au+Au collisions at y/sjfif 
= 200 GeV in the melting and default AMPT 
model with hadronic rescattering. The exper- 
imental data come from Ref. [33]. 

correlation by comparing the results of 
the default AMPT model and the melt- 
ing AMPT model. Figure 5(a) and (b) 
give three-particle correlations in selected 
(1 < A0i 2 < 5.28) regions for the melt- 
ing AMPT model and the default AMPT 
model. Note both cases are the results af- 
ter hadronic rescattering in the top 10% 
Au+Au collisions at ^snn = 200 GeV. 
Though with large statistical errors, the 
default AMPT model seems to produce a 
three-particle correlation, but the three- 
particle correlation area is considerably 
less than that from the melting AMPT 
model. It is consistent with the results 
of di-hadron correlation in our previous 
work (see Figure 5(c) and (d)) that con- 
cluded that hadronic rescattering alone 
cannot reproduce a splitting parameter 
of Mach-like structure on away side large 
enough to match the experimental mea- 
surements [21]. 



The nuclear modification factor, R cp , is also 
considered a useful probe of the energy loss 
of high pt partons in the dense medium cre- 
ated in nucleus-nucleus collisions. Figure 6 
shows the transverse momentum depen- 
dences of nuclear modification factor R cp of 
charge hadrons in the melting and default 
AMPT model with hadronic rescattering. 
The R cp is dinned by following formula: 
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E> N bin I p p T dp T dr] \ c 

~ N bin \ c x I ' 

where the Central and the Peripheral col- 
lision centralities are 0-10% and 40-60%, 
and the respective number of binary colli- 
sions N bin = 939.4 (0-10%), 93.7 (40-60%). 
The R cp in the melting AMPT model is of 
similar shape of experimental data, which 
can match experimental data well if scaled 
by a factor 1.2. However the R cp from the 
default AMPT model seems to be inde- 
pendent of pt and inconsistent with exper- 
imental data. The partonic interactions in 
the melting AMPT model appear essential 
to describe the shape of nuclear modifica- 
tion factor as a function of px in Au+Au 
collisions. Furthermore, R cp is suppressed 
more heavily in higher p? range (px > 3.5 
GeV/c) in the melting AMPT model than 
in the default AMPT model, which may 
indicate that more energies are lost into 
the medium by parton cascade mechanism 
especially for high pt particles, which is 
expected to be in favor of the formation of 
partonic Mach-like shock waves. 

V. Conclusions 

Three-particle correlations have been ex- 
tracted by using event-mixing technique 
in a multi-phase transport model with 
both partonic and hadronic interactions. 
Correlations in different azimuthal angu- 
lar regions with respect to the trigger jet 
direction, so called 'center' , 'deflected', 
'cone', 'near' and 'near-away', have been 
discussed for Au+Au collisions at ^/snn 
= 200 GeV. The AMPT results with and 
without hadronic rescattering are also 
compared. The 'center' three-particle cor- 
relation becomes more and more dominant 
with the decreasing of number of partic- 
ipants, which may reflect the centrality 
dependence of partonic density and the 
strength of partonic interactions. The den- 
sity ratio of 'deflected'/'cone' approaching 
1.0 in central collisions indicates that the 
three-particle correlation in central col- 
lisions is mainly produced by a partonic 
Mach-like shock wave mechanism, and in 
peripheral collisions deflected jet mecha- 
nism also contributes. The partonic Mach- 
like shock wave mainly originates from 



strong partonic interactions in dense par- 
tonic matter. The three-particle correla- 
tions are also sensitive to hadronic rescat- 
terings, therefore the effect of hadronic 
rescattering may need to be considered in 
quantitative studies. The default AMPT 
model, where only the hadronic rescat- 
tering mechanism plays a dominant role, 
produces a three-particle correlation area 
much smaller than the melting AMPT 
model which includes both parton cas- 
cade and hadron rescattering mechanisms. 
Our AMPT calculation of three-particle 
correlations re-affirms our previous con- 
clusion from di-hadron correlation studies 
that hadronic rescattering alone cannot 
produce an amplitude of Mach-like cone 
on away side large enough to match the 
experimental data. Parton cascade mecha- 
nism is essential and important in order to 
describe the amplitude of observed experi- 
mental Mach-like structure. 
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